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i Abstract 

We generalise Livens theorem, showing that Hamiltonian equation on the vector bundle 
E* — > M, dual to a general algebroid E — > M, can be derived by means of a variational 
principle. The framework can be used to describe Hamiltonian dynamics associated with 
(both vaconomic and nonholonomic) constraints in the bundle E. 
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1 Introduction 



■ Lagrangian mechanics on general algebroids The importance of a Lie algebroid structure 

in the formulation of Geometric Mechanics is now commonly accepted. The basic reason is that 
| Lie algebroids cover, on the one hand, the structure of the tangent bundle of a manifold and, 

on the other hand, geometric structures which arise from a reduction by a symmetry group such 
as a Lie algebra of a Lie group or, more generally, a space TP/G of G-invariant vectors on a 
principal bundle G — > P — > M. Consequently, the formulation of Geometric Mechanics on Lie 
algebroids allows one to treat many natural situations in frames of a single universal formalism 

Such a formulation was proposed by many authors [2j [TTJ [9j [TOJ EI2 E] • In this paper we 
will base on an approach by Grabowska, Grabowski and Urbanski [H [5] , in which the dynamics 
associated with an algebroid E is derrived by means of a certain double vector bundle morphism, 
naturally generalising the celbrated Tulczyjew triple |17j . These authors were also the first 
to realise that neither the Jacobi identity, nor the skew-symmetry of the bracket is necessary 
to generate the dynamics. Consequently, their framework can be formulated also for general 
algebroids (cf. [6]). Such a generalization proved to be useful, for instance, in the context of 
nonholonomic constraints [8]. 

A natural problem is to incorporate the constraints into the framework of geometric mechanics 
on algebroids. This can be relatively easily done for a Lagrangian formalism. Namely as shown in 
[5] the generalised Euler-Lagrange equations can be derived by means of variational principle for 
E- valued curves. Variational calculus is here understood as the study of the differential dSii'j) 
of an action ^ 

5 L ( 7 ) := / 1 L(l(t))dt, 
Jt 
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where 7 : [to>£i] — > E is an .E-valued curve. For a variation £7 G T 7 £" we simply have 

< ( 5 7 ,d < S L ( 7 )> := r<5 7 (t),dL( 7 (t))>dt. 

Our attention is restricted to certain curves 7, called admissible trajectories and certain variations 
^7, called admissible variations, which can be naturally constructed from an algebroid structure 
on E (see section [3] for details). Constraints can be now understood as restrictions in the set of 
admissible variations and/or trajectories. Constrained Lagrangian dynamics on algebroids can 
be now derived from a variational principle for these restricted sets [5]. 

For the constrained version of the Hamiltonian dynamics there is so far no such an easy 
approach, even if an algebroid is simply a tangent bundle. Typically the constrained Hamiltonian 
equations are derived from the constrained Euler-Lagrange equations via coordinate calculations 
under an additional assumption of hyperregularity of the Lagrangian. This has been done by 
|15] for the case of vaconomic constraints and by [1] [16] for nonholonomic constraints. 

Livens theorem In this paper we propose another method of generating the constrained 
Hamiltonian dynamics on general algebroids. Our motivation comes from a classical Livens 
theorem stating that the unconstrained Hamiltonian dynamics on T*M has a variational for- 
mulation. Namely, the Hamiltonian trajectories can be described as the critical points of the 
action ^ 

S H (x(-),p(-)) := [ 1 (x(t),p(t)) - H(x(t),p(t))dt, 
Jt 

with respect to variations (5x(-),5p(-)) induced by homotopies in T*M. (We use standard no- 
tation (x a ) for coordinates on M, (x a ,x b ) and (x a ,pb) for natural coordinates induced on TM 
and T*M). 

Our results If we reformulate the above results in terms of the TM (Bm T*M-trajectories 
(x(-),y(-),p(-)), namely as a study of the action 

S H (x(-),y(-),p(-)) := f 1 (y(t),p(t)) - H(x(t),p(t))dt, 

under an additional admissibility condition y(t) = x(t), Livens theorem can be easily generalized 
to the context of general algebroids as is done in theorem 15.11 To be more precise, given a 
Hamiltonian function H : E* — > R, we can define a function Ch : E (Bm E* — > R by the formula 

C H (Y) := (p E (Y),p E *(Y)) -H(p E *(Y)), 

where p E : E(BmE* — > E and p E * : E®mE* — > E* are natural projections. Now for a trajectory 
r : [to, t\] — > E (Bm E* we can define the action 

S£ H (r) := t C H (T(t))dt. 
J t 

We restrict our attention to trajectories T such that 7 := p E oT : [icb^i] - > E is an admissible 
trajectory in E (in the sense of |5]) and variations 5T : [to, t±] — ) Tr(E (Bm E*) ~ T Pe ^E (Btm 
T p »(T)E* such that p TS o5T : [to,ti] — >■ T-yE is an admissible variation in TE (in the sense 
of [5]). Extremals T of the action Sc H (-) under these restrictions give Hamiltonian dynamic 
on E* defined by a vector field Xh := tdH^E, where n^; is a contravariant 2-tensor on E* 
equivalent to the presence of an algebroid structure on E. Precisely, the .E*-part p E * oT is a 
Hamiltonian trajectory, and the E-pait p E oT = T*tt o dH can be understood as a momenta - 
velocity correspondence (Legendre map). 
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Constrained Hamiltonian dynamics As we mentioned above, the derivation of the con- 
strained Lagrangian dynamics on an algebroid E was possible by restricting the sets of admissi- 
ble variations and trajectories. On the other hand, the (unconstrained) Hamiltonian dynamics 
is obtained via a variational approach which uses the unconstrained admissible variations and 
trajectories on E. A natural idea is hence to use the variational principle for the Hamiltonian 
dynamics for sets of admissible variations and trajectories restricted by the constraints. In con- 
sequence, we obtain a constrained version of the Hamiltonian dynamics on general algebroids. 
The results agree with the equations derived by |15| [TJ [16] , although in contrast to these papers 
we start from an arbitrary Hamiltonian function, not assuming hyperregularity. 

Organization of the paper After introducing general algebroids in section [2j we sketch an 
abstract and quite a general approach to variational calculus (section |3|) . In section [4] we show 
its particular realization in the case of Lagrangian dynamics on a general algebroids (following 
the original paper [5j). In the same section we also define Hamiltonian dynamics on algebroids. 
Section[5]is devoted to the formulation and proof of theorem 15 . 1 1 which generalizes Livens theorem. 
Finally in section we calculate the constrained Hamiltonian dynamics in the vaconomic and 
nonholonomic case and compare our approach with other authors. 

2 General algebroids 

Local coordinates, notation Let M be a smooth manifold and let (x a ), a = l,...,n, be 
a coordinate system in M. We denote by tm '■ TM — > M the tangent vector bundle and by 
ttm '■ T*M — )• M the cotangent vector bundle. We have the induced (adapted) coordinate systems 
(x a ,x b ) in TM and (x a ,pb) in T*M. More generally, let r : E — > M be a vector bundle and let 
7r : E* — > M be the dual bundle. Let (ei, . . . , e m ) be a basis of local sections of r : E — >■ M 
and let (e\, . . . ,e™) be the dual basis of local sections of tt : E* — >■ M. We have the induced 
coordinate systems: (x a ,y l ),y r = i e i in E, and (x a , £j) , £j = i e . in E*. 
Consequently, we have natural local coordinates 

(x a ,y\x b ,f ) in TE, (*°, m TE*, 

(x^y^p^TTj) in T*E, (x», in T*£*. 

For a curve j(t) ~ x a (t) in M we will denote its tangent prolongation by T(j(t)) ~ (x a (t), x b (t)) 
in TM. 

General algebroids An algebroid structure on a vector bundle r : E — > M is given by a 
bilinear bracket [•, •] on the space Sec(E) of (local) sections of r, together with a pair of vector 
bundle morphisms p,a : E — > TM (left and right anchor maps) such that the compatibility 
condition 

[g-X,f-Y] = gf[X, Y] + g ■ p(X)(f)Y - f ■ a{Y)(g)X (2.1) 

is satisfied for every X,Y G Sec(£) and f,g G C°°(M). 

In literature one considers often certain subclasses of algebroids. Under additional assumption 
that the bracket is skew-symmetric (hence p = a) we speak of skew- algebroids. A skew-algebroid 
whose anchor map is an algebroid morphism, i.e. 

p([X,Y]) = [p(X),p(Y)] TM (2.2) 

is called almost-Lie algebroids. Almost-Lie algebroids with the bracket satisfying the Jacobi 
identity (in other words: if the pair (Sec(2£),[-, •]) is a Lie algebra) are called Lie algebroids. 
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In the context of mechanics it is convenient to think about an algebroid over M as a gener- 
alisation of the tangent bundle of M. An element a € E has an interpretation of a generalized 
velocity with actual velocity v 6 TM obtained by applying the left anchor map v = p(a). 

In local coordinates (x a ,y l ), as introduced at the beginning of this section, the structure of 
an algebroid on E can be described in terms of function pf(x), crf(x) and c l - k {x) on M given by 

p(ei) = pf(x)d x a, cr(ej) = af(x)d x a and [e»,ej] = c^(x)e k . 

Examples 

The tangent bundle TM. Every tangent bundle tm '■ TM — > M of a manifold M posses 
a natural structure of a Lie algebroid. We simply take the anchor map to be idxM : TM — > TM 
and the standard Lie bracket of vector fields on M as the algebroid bracket. 

The Lie algebra q of a Lie group G. A Lie algebra q understood as a vector bundle 
over a single point also posses the structure of a Lie algebroid. The anchor map is in this case 
trivial, whereas the algebroid bracket is simply the Lie bracket [■, -] . 

The Atiyah algebroid TP/G. An example generalizing the above two arises from the 
study of G-invariant vector fields on the total space of a principal G-bundle G->P4 M. Every 
such field can be canonically identified with a section of the quotient bundle E = TP/G — > M. 
Since a Lie bracket of two G-invariant vector fields is again G-invariant, Sec(E) is equipped with 
a natural Lie algebra structure (Sec(E), [•, -]e)- This, together with the natural anchor map 
p : E = TP/G —> TM defined as p([X]) = ir*X, forms a Lie algebroid. It is known as the Atiyah 
algebroid. 

The above are examples of Lie algebroids. Natural examples of skew-algebroids appeared 
first in the context of nonholonomic constraints (cf. for instance [5]), while the definition of a 
general algeboid was introduced in |6], where some natural examples are given. 

The structure of an algebroid has several equivalent descriptions. We will give a few of them 
which are particularly useful in variational calculus - compare [5J S] . 

Algebroids as linear (2,0)-tensors on E* The algebroid structure on E corresponds to a 
presence of a linear (2,0)-tensor field II on E*. In local coordinates 

n = 8) % + p\{x)d^ 8) 8 xb - a b i (x)d Xb %. 

The linearity of II means that the natural morphism tn : T*E* — > TE* induced by the contraction 
is a morphism of double vector bundles. 

T*E* ^TE* (2.3) 




For Lie algebroids II is a linear Poisson structure on E*. It is well recognised in two standard 
examples: if E = TM, II is a canonical Poisson structure on T*M, whereas if E = g it is a 
Lie- Poisson structure on q*. 
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Algebroids as double-vector-bundle morphisms Composing double vector bundle mor- 
phism ljj with a canonical double vector bundle isomorphism (and antisymplectomorphism) 
R T : T*E ~ T*^* one obtains double vector bundle morphism ee '■ T*E — > TE*. The presence 
of such a morphism provides an equivalent characterisation of the algebroid structure on E. 

T*E >■ TE* (2.4) 




Since in local coordinates R T is given by 

R T (x a ,y\p b ,TT j ) = {x a 1 TT i ,-Pb,y j ), 
the morphism ee is of the form 

e E {x\y\p h ^) = (x a ,^,p b k ^)y k ,4^)y% + ^)Pa)- (2.5) 

Relation k Since the dual bundles of tte ■ T*E — >■ E and Tit : TE* — > TM are, respectively, 
te '■ TE — > E and Tr : TE — > TM, the dual to eg is a relation ke ■ TE — > TE. It is a uniquely 
defined smooth submanifold k in TE x TE consisting of pairs (v, v') such that p(te(v')) = Tt(v) 
and 

(v,e e {v*)) Tt = (v',v*) te 

for any v* € T*_ e ^ v ,^E, where -} Tt is the canonical pairing between TE and TE*, and <^-, 

is the canonical pairing between TE and T*E. We will write k : v — > v' instead of (v, v') € k. 

This relation can be put into the following diagram of "double vector bundle relations" 

TE •* TE (2.6) 




In local coordinates it reads 

« : (>, Y\ p b k (x)y k , yi) - (x a ,y\ a b k {x)Y k , ^ + ^(x)/r') . (2.7) 

A canonical example of a mapping in the case of E = TM is given by £tm = ~ ^ ne 
inverse to the Tulczyjew isomorphism olm ■ TT*M —> T*TM [17J. The dual relation is in this 
case the well-known "canonical flip" km '■ TTM — > TTM. Since au is an isomorphism, km 
is a true map, in fact - an isomorphism of the corresponding two vector bundle structures as 
well. In local coordinates (x a ,x b ,5x c ,5x d ) on TTM, km maps an element (x a ,x b ,5x c ,5x d ) to 
(x a , 5x b , x c , 5x d ). 
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3 An abstract approach to variational calculus 



Variational problems A standard data in a variational problem consists of a vector bundle 
r : F —> M (perhaps with additional structure) and a function L : F — > R (Lagrangian). Given 
any smooth path 7 : I = [to, ti] — > F we can define the action 

S L ( 7 ) :=^L( 7 (t))dt. 

In typical situations we do not want to consider Sl for all 7 6 C°°(I,F), but instead we restrict 
our attention to a certain subset T C C°°(I, F). Curves in T will be called admissible trajectories. 
As an example consider the tangent bundle tm '■ TM — > M. Typically we are interested only in 
curves 7 : I — > TM which are tangent prolongations of curves in M, i.e. j(t) = T (tm 7 (t))- 

A variation (or a virtual displacement) of 7 € T is a curve £7 : I — > Ti* 1 which projects to 
7 under tf '■ TF -4 -F. The set of all variations will be denoted by C. We will use a symbol 
C 7 for the set of all variations which project to 7. Among elements of C we can distinguish a 
class C° of variations with vanishing end-points. Denote by C® the common part C°nC 7 . A earn, 
sometimes we will like to restrict our attention to admissible variations only that is variations 
from a certain subset W C C. We will use a natural notation W° := W n C°, W 7 := W n C 7 and 
W° := wnc°. 

For a given 7 G T and £7 S W 7 we can consider the variation of the action Sl at 7, along 

£7 

(<57,d5 L (7)>:=^<57(t),dL( 7 (t))>dt. 

We would like to understand a variational problem as the study of the critical points of the action 
Sl restricted to curves in T and variations in W. The variational problem can be hence defined 
as a quadruple 

V :={F,L,T,W). (3.1) 

Of course, in typical situations the sets T and W are not arbitrary, but are somehow related to 
the geometry of the bundle F. 

A curve 7 € T will be called an extremal of V iff 

(5 7 , dS L ( 7 )} = for every <5 7 G W°, 

i.e. 7 is a critical point of the action Sl relative to admissible variations with vanishing end- 
points. The set of all extremals of V will be denoted by T-p C T ■ 

Constraints The above definition of the variational problem is very general and particularly 
useful in the context of constrains. We can define V - a constrained variational problem by 
restricting the set of admissible trajectories and/or admissible variations, i.e. V = (F, L,T , W ), 
where "Y C T and W C W. In this way we follow the ideas of Tulczyjew [18], namely that 
constraints should be imposed rather on virtual displacements than on configurations. Usually 
these restrictions are somehow related to additional geometric structures in the bundle F . We will 
see concrete examples while studying the nonholonomic and vaconomic constraints associated to 
a subbundle D of an algebroid E. 

4 Lagrangian and Hamiltonian dynamics on algebroids 

In section [3] we defined an abstract variational problem as a quadruple, consisting of a vector 
bundle, a Lagrangian function and the sets of admissible trajectories and admissible variations. 

Now we will show that if the considered vector bundle is equipped in a general algebroid 
structure, then the sets of admissible trajectories and variations can be defined in a natural way. 
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The algebroid variational problem T*e Consider a vector bundle r : E — > M with a 
structure of a general algebroid. Among the curves 7 : [toj^l] E we distinguish admissible 
paths (or E-paths shortly), such that the tangent prolongation T(x(t)) of their base curves 
x{t) := r o 7(f) coincide with their left anchor 

T(x(t)) = p( 7 (t))- (4.1) 

The set of all such curves will be denoted by Te- 

Let now 7(4) € E be an E'-path over x{t) G M and let 6(f) £ £7 be any other path covering 
x(t). The tangent prolongation T(b(t)) is a vector field in TE along b(t). There exists an unique 
vector field <5&7(i) in TE along 7(f) (a variation of "f(t)) which is K-related to T(b(t)) (cf. [5] for 
details). In local coordinates in TE, if -y(t) ~ (x a (t) , y { (t)) and b(t) ~ (x a (f), &*(£)), then 

8 bl (t) ~ (^(^.^(^.^(^(t),^^) +4(a:)y i (t)y'(t)) . (4.2) 

The set of variations of the above form will be denoted by We- We may say that the variation 
<^> 7 (i) is generated by 6(f) - the infinitesimal variation. The subset Wg C We will consist of 
variations generated by b(t) 's vanishing at the end-points. 

Now, for a given Lagrangian function L : £7 — > R we can consider an algebroid variational 
problem 

Ve:=(E,L,T e ,W e ). 

Admissible paths and admissible variations have a nice geometric interpretation in the context 
of a Lie grupoid - Lie algebroid reduction. Namely, admissible paths are paths in a grupoid 
reduced to the algebroid and admissible variations are infinitesimal homotopies in a grupoid 
reduced to the algebroid. For a detailed discussion the reader should confront [3J |7]. Let us 
mention only that the standard variational problems are of the above type. 

Classical variational problem on TM Consider the standard (unconstrained, with fixed 
end-points) variational problem of finding the extremal x(t) G M of the action 

h 

L(x(t),x(t))dt, 

where L : TM — > R. This problem is, in fact, a variational problem on a vector bundle tm '■ 
TM — > M with admissible trajectories being just the tangent prolongations of the base curves 

7(t) ~ (x(t),v(t)) G TM, where v(t) = x(t). 

Admissible variations are defined by means of base homotopies 

5j(t) ~ (5x(t),6v(t)) = (d s x(t,s)\ s=0 ,d s x(t,s)\ s=0 ) , 

where x(t,s) is the homotopy in M and x(t, 0) = x(t). In fact every such variation is generated 
by a curve b(t) = d s x(t, s)\ s= q G T x mM. Indeed, 

(5x,5v(t)) = (b(t),K M (b(t))) , 

where km '■ TTM — > TTM is the canonical flip which maps d s (dtx(t, s)) into dt (d s x(t, s)). The 
presence of km on TM is often not recognised, being hidden under the obvious symmetry of 
second derivatives. In local coordinates (x a ,x b ) on TM, if b(t) ~ (x a (t) , b c (t)) , then 

5x a (t) = b a (t) and 5x c = b c (t). 

We recognise equation (|4.ip and (|4.2p for the standard algebroid structure on TM (in coordinates 
(x a ,x b ) on TM we have clearly p^(x) = cr"(x) = 5" and c a de {x) = 0). 
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Euler-Lagrange equations For the algebroid variational problem Ve = (E, L,7~e,We) one 
can derive an analog of Euler-Lagrange equations in a simple manner [5J [5]. Take 7 6 Te and 
(^7 G (W£;) 7 . Let us calculate 



(<5 6 7,dS L ( 7 )) 



fi 



(* 6 7(t),dL(7(t))>d*. 



(4.3) 



In local coordinates (x a ,y t ) on £" denote 7(t) ~ (x a (t), y l (t)) and 6(i) ~ (x a (t) , b l (t)) . Now 



* 67 (t),dL( 7 (i))) = a?(xW(t) — ( 7 (t)) + (6* + C ;. fc (*y (t)6 fc (i)j ^( 7 (i)) 



ax 



i9f r)f rl r)^ 

°f W^(7(0) + 4(x)y*(t)^( 7 W) " ^^(7(*)) 



<9y A 



dt dy % 



d 

+ d7 



6^L l ( 7 (t)) + ^[6 l (t)(A L ) J ( 7 (t))], 



where (Al)«(x, y) := §±(x,y) and 5L;(a;,y) := of(x)^(x, y) + e£ {x)y^{x, y) - &§jp(x, y). 
Consequently, 



<<5 6 7(i),d5 i ( 7 (i))) = (b{t)AMt)))\% + / (6(t),5L( 7 (t)))dt. 



to 



Now, j(t) is an extremal of Ve, if the above expression should be zero for all 6's vanishing at 
the end-points. This gives <5L( 7 (i)) = 0. In local coordinates 7 (i) ~ (x a (t), y l (t)) that reads 



H r)T r)T r)T 

- — ^{t),y{t))=at{x{t))—^^ (4.4) 



which, together with an admissibility condition for 7 (/j) 

x a (t)=p?(x(t))y i (t), 



(4.5) 



generalises the Euler-Lagrange equations. 

The above equations can be also obtained in a purely geometrical way using the diagram 
(|2.4p (confront [U [5]). Given a Lagrangian function L : E — > M. one can construct two maps 
Al = te* dL : E E* and = ee° dL : E — > TE* as shown on a diagram 



T*E ■ 

dL 

E- 



TL* 



In local coordinates (cf. (j2.5j) ) 

Ai(x,y) ~ ^x a , ^( x ,y) 
dL 



dL 



dL 

A L (x,y) ~ ( x a ,^ I (x,y),p b l (x)y\c^ j (x)y t -^(x,y) +^{ x )q^ (*,v) 



Now, Euler-Lagrange equations f|4.4[> and (|4.5p are simply equivalent to A£( 7 (f)) being the 
tangent prolongation of Ai,( 7 (i)) 



T(A£( 7 (*))) = A L ( 7 (t))- 
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Vaconomic constraints Let now E be an algebroid and D C E be a vector subbundle (over 
a full base M - compare remark [673]) . To define a vaconomically constrained variational problem 
on an algebroid E, associated with the bundle D, we restrict our attention to admissible paths 
in D and variations tangent to D. That is, 

Td,v<ic := {7 e Te : t(*) g -° for ever y t g [*o 3 *i]} and 
WD,vac ■= {Sb7 G W E : tf 6 7(t) G T 7(t) D for every t € [*o,*l]}- 

In the standard case E = TM, this variations correspond simply to infinitesimal homotopies of 
paths tangent to the distribution D. Note that if ^7 € Wd,vo,c, we have a priori no control of 
the infinitesimal variation b(i). Since the variations are tangent to D, the vaconomic dynamics 
is determined by the restriction of L to D. 



Nonholonomic constraints With the same subbundle D C E one can associate also a non- 
holonomically constrained variational problem. In this situation we restrict ourselves again to 
admissible trajectories in D 

,vaci 

but input different restrictions on the set of admissible variations. Namely, according to 
d'Alembert principle, we consider only these variations 5y) £ We which are generated by in- 
finitesimal variations b(t) with values in D 

yV D ,nh = {hi G Wi? : b(t) € D for every t € [t ,h]}. 

In this case we have a full control of the set of infinitesimal variations b(t). Note, however, that 
the variations 5&7(£) are, in general, no longer tangent to D. For this reason the knowledge of 
L\d is not sufficient to study nonhonomically constrained dynamics. 

In cases of both, nonholonomic and vaconomic constraints associated with D, the study of 
action (|4.3|) with admissible variations restricted to yVo,nh an d V^D,vac, respectively, allows one 
to derive the algebroid version of the constrained Euler-Lagrange equations on E (confront |5]) 
which generalizes the standard constrained dynamics on TM. 



Hamiltonian dynamics As we have already mentioned, an algebroid structure on E is equiv- 
alent to a presence of a (2,0)-linear tensor field He on the dual bundle E* (linear Poisson bivector 
field for Lie algebroids). Such a structure allows one to define a natural analog of the Hamiltonian 
dynamics on E* . 

Namely, for a given function H : E* — > R we can define a Hamiltonian vector field Xjj in E* 
by a natural formula Xjj := tdH^- Explicitly, the associated dynamics of Xjj is governed by the 
following equations 



dH 



dH 

' dx a 



(4.6) 



As in the case of the Euler-Lagrange equations it is convenient to describe the construction 
of Xjj by means of a diagram (|2.3p 



T* e* 

dH 

E* - 



■ TE* 



id 



■E* 
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In the next chapter we will derive Hamilton equations (|4.6p by means of a variational principle 
in the sense introduced in previous sections. What is more, the general approach to constraints 
presented before allows one to derive Hamiltonian equations on F* associated with constraints 
in the configuration bundle F. 

5 Variational description of Hamiltonian dynamics 

Geometric setting Let us introduce a few definitions. We will develop a variational calculus 
on a bundle F := E ®m E* M. The tangent bundle TF of F can be canonically identified 
with TF = T(F M F*) w TE ® TM TF*. Both F and TF admit natural projections which will 
be denoted by p E : F — > E, p E * : F — > E* , p TS : TF — > TF and p^E* : TF —> TF* respectively. 
Finally, we have natural coordinates (x a ,y l ,^j) on F and (x a , y l , £j , x b , y k , £i) on TF. 

For F as above we can construct sets of admissible trajectories Tf arid admissible variations 
Wf as follows 

Tf := {r : [t , h] -> F : p E oF € T E } , 
W F := {ST : [t , h\ -> TF : p TE o5T e W E } ■ 

That is, admissibility is determined by the F or TF part of a path or a variation only. We will 
distinguish the class W F consisting of those variations SF which project to W E under p^ E . 

For a given function H : F* — » R (Hamiltonian) we may define a new (Lagrangian) function 
Ch '■ F — >• R by the formula 

C H (X) := < Ve X, p E * X) - H( Ve , X). 
It defines an action on the space of admissible trajectories in F 

Sc H (T) = f 1 (v E or(t),p E * oT(t)) - H( VE « oT{t))dt. 
J t 

Now we have all ingredients (namely F, Tf, Wf and Ch) necessary to develop variational 
calculus on F in the sense of section [3] 

Theorem 5.1. The trajectories £(t) € F* of the Hamiltonian vector field Xh are precisely the 
F* -projections of the extremals V of the action Sc H { ) with respect to variations in W E . 

For critical Y{t) which projects to £(t) € E* , the E-projection is simply p^(r(t)) = Ttt o 
dH(Z(t))€E. 

Proof. For simplicity we may assume that the extremal T(t) ~ (x a (t) , y l (t) , £j (t)) lies in a single 
coordinate chart of F ®m F*. 

Every admissible variation 6T G Wf along T has the form 

ST(t) = a?(x)f(t)d xa + (f k + c^xWWfit)) d yk + h k (t)dt„. 

Now, 

<5r,d5 £H > = ( VTE o5r(t),p TE * o5T(t)) TT -(p TE , o5T(t),dH(p E , oT(t)))dt = 

J t 

rti flj-f f)ff 

= J (f k + c%{x)y\t)f {t))Ut) + y k (t)h k {t) - of ( x )/*(t)— - h k (t)—dt. 
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Integrating by parts we will get 



[ST,dS CH ) = fit 



If ST G Wp, then vanishes at to and t\, hence 

(ST, dS CH ) = jT 1 / J ' (-6 " ^(*)f| + &4 + h k (t) (y k - ||) dt. (5.1) 

The last expression vanishes for every such 5T if and only if 

dH 

i j = ^ J (x)y i -^(x) — (x,0 

and 

i, . ... 

Consequently, since by assumptions p E oT(t) ~ (x a (t) , y l (t)) is admissible, that is x a = pf(x)y l , 
we get 

f) J-f 

in perfect agreement with ([4.6]) . In other words, p E * oT(t) ~ (x a (t),^j(t)) is a Hamiltonian 
trajectory of Xh- What is more, p s » or (f) ~ (x a {t), §^(x(i), £(£))), which is precisely the 
coordinate form of T*ir o dH(£(t)). □ 



6 Hamiltonian dynamics with constraints 

As we have seen, Hamiltonian dynamics can be described by means of a variational principle on 
F = E® M E*, provided that we impose simple condition on admissible variations p^E 
and trajectories p E oT E Te- It is a natural idea to define constrained admissible variations and 
trajectories, associated with a subbundle D C E, in a similar manner; 

Tf,d,vo.c := {r G Tf : p E or 6 Td,vo,c}, 

TF,D,nh ■= {r £ Tf ■ Y>E °r £ TD,nh}, 
W F ,D,vac ■= {ST e W F : PTE °ST € }V D , TC }, 

w F ,D,nh := {^r g w F ■ pte e w D)nh }, 

and to derive a constrained version of the Hamiltonian dynamics using the variational approach. 
It seems also possible to consider constraints in the phase part E* of F. However, since E* does 
not posses a priori any natural structure to generate constrained variations (such as an algebroid 
structure), these could be only configuration constraints. 
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Nonholonomic constraints In this case we can calculate the associated equations of motion 
in a manner similar to the proof of theorem 15.11 Indeed, we can derive formula (|5.ip again, 
yet now, since 5T £ Wp,D,nh the curve f l (t) belongs to D, whereas /ifc(i) can be arbitrary. 
Consequently, we get that 

dH 

annihilates D. In other words, if D is given locally via implicit equations <& s (x,y) = 0, there 
exists functions fJ, s (t) such that 

i 3 = ^(x)y* - ^(x) — (x, + »s(t)-Q- [x, — 

In addition, we have the condition y k = §jr(%,fi) an d the fact that the E'-projection of the 
extremal p E oT ~ (x a ,y*) = (x a , jjjr) should belong to Tb.nh- As a consequence we get the 
following set of differential equations 



• ^ kt ,dH ( ar .dH ,.d<S> s ( dH 

f) H 

x a = P a k (x) — {x,t), (6.2) 
Oc,k 

d H 

^(*,_) = 0. (6.3) 

The same equations (yet for the standard Lie algebroid structure on TM) have been ob- 
tained for instance in [16j. In the derivations the authors assumed that the Hamiltonian comes 
from a regular Lagrangian via the Legendre map. In our approach no such an assumption is 
necessary - we may simply start from an arbitrary Hamiltonian. Equations (|6.ip - (|6.3p . again for 
Hamiltonian given via a Legendre map of a regular Lagrangian, can be also obtained from the 
nonholonomically constrained Euler-Lagrange equations on algebroids - compare |5] (equations 
(5.10) and (5.11)). 

Geometric approach to nonholonomic constraints Equations fj6 . 1 1) — fj<3. 3|) can be described 
in a purely geometric way. The construction is based on the diagram (|2,3p . Observe first that 
the inclusion i : D — > E gives rise to a dual map i* : E* — > D* and consequently to a tangent 
map Ti* : TE* — > TD* over i* . Composing (|2.3p with the later we will obtain the following 
diagram 

T?'* 

T * E * TE* 1 TD* , 



dii 



*2 - 



TD* 



E* -~"~ id ^ E* D 



where X§ := Ti* o X H is a map from E* to TD*. 

Theorem 6.1. The curve j(t) ~ (x a (t), £i(t)) £ E* satisfies the nonholonomic Hamilton equa- 
tions (|6.ip - (|6.2p if and only if XS^y(t)) G TD* is a tangent prolongation of a curve i*7(t) € D* . 
Condition (|6.3p means simply that T*7r(dH) belongs to D. 

Proof. For X € TE* consider a vector 

6H(X) := X h {te*{X)) -Xe TE*. 
In local coordinates if X ~ (x a ,^i,x b ,S,j) then 

(dH dH dH \ 
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Now X = T(7(i)) satisfies (|6.2p if and only if 5H(X) is a vertical vector, whereas (|6.ip means that 
its vertical part V(5H(X)) =: SH(X) ~ (>, £ fcC f .(^f^, - <7*(x)J^(s,0 - &) belongs to 

the annihilator D° C E*. Observe that 6H(X) € D° if and only if i*5H(X) = 0. Finally note 
that since i* : E* — > D* is the identity on the base, a vector Y € TE 1 * is vertical if and only if 
Ti*Y is vertical. Moreover, for a vertical vector Y £ FS* we have i*V(Y) = V(Ti*Y). 

Collecting all the above fact we conclude that (j6.ip and (|6.2p is equivalent to 5H(T(j(t))) £ 
FE 1 * and i*V(£ff(T(7(t)))) = which is, in turn, equivalent to Ti* (?iJ(T(7(i)))) = 0. This, 

by definition of <5iJ and Afjy, means that 

XgW)) = Ti*X H W)) = T(i* 7 (t))- 

□ 

Nonholonomic reduction Considerations from the last paragraph suggest that the dual 
bundle D* plays the role of a phase space for the nonholonomic dynamics. This in turn leads to 
a natural question is it possible to formulate the dynamics purely in therms of the bundle D*. 
Such a procedure is known as a nonholonomic reduction. 

This problem for D C TM was considered by van der Shaft and Maschke [16J. The authors 
describe the nonholonomically constrained Hamiltonian dynamics on T*M (for regular Hamil- 
tonians) as an unconstrained Hamiltonian dynamics on D* obtained using a "Poisson" bracket 
on D* which needs not to satisfy the Jacobi identity. Using our terminology the structure they 
consider is simply a bi-tensor field on D* which corresponds to a general algebroid structure on 
D. 

The above results seems very interesting as they allow to cover the constrained and un- 
constrained dynamics within a universal algebroid formalism. Unfortunately the nonholonomic 
reduction of |16| has no geometric meaning. Strictly speaking the authors are unaware of using 
the decomposition TM = D © D which depends on the choice of local coordinates. 

In some cases, however, where a natural decomposition E = D © D can be given, the 
nonholonomic reduction can be carried out. Below we will state a result of this kind, which can 
be understood as a Hamiltonian analog of the Lagrangian nonholonomic reduction described in 
[8]. We do not give a proof (which is quite simple) since this paragraph was intended to be 
informative only. 

Theorem 6.2. Consider a nonholonomically constrained algebroid variational problem associated 
with a subbundle D C E of an algebroid (E, [•, •],/>, a) for a Hamiltonian H : E* — > R. Assume 
that H is of mechanical type, i.e. H(x,£) = hg~ l {£,£) + V{x), where g is a non- degenerate 
symmetric bilinear form on E and V : M — )■ R is a smooth base function. 

Let E = D © D 1 - and E* = D* © (D -1- )* be natural decompositions and P : E — > D a natural 
projection induced by g. 

The nonholonomically constrained Hamiltonian trajectories on E* are the same as uncon- 
strained Hamiltonian trajectories on D* derived for a Hamiltonian i?|.o*©{0} an d an algebroid 
structure (D, [•, -]d, pu, gd) given by the bracket [-,-]d '■= •P[v]> an ^ anchors := p\o and 
&d '■= o~\d- 

Vaconomic constraints The derivation of the constrained equations of motion is much more 
delicate in this case. Repeating the derivations from the proof of theorem 15. II we can again obtain 
(|5.ip . We can easily conclude that y k = ^-(x,^), since hk{t) can be arbitrary. Now, however, 
we have no clear condition for f k (t) apart from p TE oST € TD, so it is difficult to derive the 
second equation. 
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Instead, we can give a simple sufficient condition for T being the extremal in this case. 
Observe that if <I> : IR x .E — >• M is a function vanishing on D, then, since p TE o<5r is tangent 
to D, we have {5T(t),d$(t,p E oT(t))\ = 0. Now, if for any such <I> the curve T is an extremal 
of the unconstrained variational problem on F for the new Lagrangian function £h(X) := 
C H (X) + $(t,p E X), then for ST € W° F ^ vac we have 

= (ST(t),d£ H (T)(t)) = (5T(t),d£ H (T(t)) + d$(t,p E oT(t))) = (8T(t),dC H (T(t))). 

In consequence, T is an extremal of the vaconomically constrained problem. If D is given by 
the set of implicit equations <& s (x,y) = 0, we may take $>(t,x,y) of the form fj, s (t)Q s (x,y). The 
equations of motion in this case read 



d/ , .d$"\ (, . ,d$ s \ k , .8H ar ,(dB . 

F) ff 

i a = p a k (x) — (x,0, (6.5) 

with the additional set of conditions: 

<S>»(x,—) = 0. (6.6) 

To link the above equations to known results, note that the vaconomicaly constrained Euler- 
Lagrange equations on general algebroids have been derived in j5] (equations (5.3) and (5.4)). In 
case that the Hamiltonian H comes from a regular Lagrangian L via a Legendre map, equations 
(|6.4p - (|6.6p can be eaisily derived from these equations. 



Remark 6.3. In this paper we introduced the constraints (both nonholonomic and vaconomic) 
by means of a subbundle D C E over a full base M. However, most of the results would 
remain valid also if D was a subbundle over a submanifold Mo C M if we assumed admissibility 
conditions p(D) C TM^ and cr(D) C TMjj. We decided to work in a less general situation for 
notation simplicity. 
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